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We explain the findings by Di Leonardo et al. [Phys. Rev. Lett. 84, 6054 (2000)] that the 
effective temperature of a Lennard- Jones glass depends only on the final value of the density in the 
volume and/or temperature jump that produces the glass phase. This is not only a property of the 
Lennard- Jones liquid, but a feature of all strongly correlating liquids. For such liquids data from 
a single quench simulation provides enough information to predict the effective temperature of any 
glass produced by jumping from an equilibrium state. This prediction is validated by simulations 
of the Kob-Andersen binary Lennard- Jones liquid and shown not to apply for the non-strongly 
correlating monatomic Lennard- Jones Gaussian liquid. 



Condensed matter is frequently found in out-of- 
equilibrium states. For example, for systems like super- 
cooled liquids, dense colloids, spin systems, etc., the (off- 
equilibrium) glass state occurs naturally after cooling or 
compression from a state of thermal equilibrium. An ef- 
fective temperature describes the non-equilibrium prop- 
erties of a glass, and the possibility of connecting the 
effective temperature with the observed violation of the 
fluctuation- dissipation theorem 1 (FDT) has opened new 
ways of inquirjJ^lED In 2000 Di Leonardo et al. 5 stud- 
ied the off-equilibrium dynamics of the single-component 
Lennard- Jones (LJ) liquid (with a small many-body term 
added to the potential to prevent crystallization). This 
system was subjected to sudden temperature decreases at 
constant density (quenches) as well as to sudden density 
increases at constant temperature (crunches). From the 
violation of the FDT the effective temperature was de- 
termined. Surprisingly, it was observed that the effective 
temperature T e ff is independent of the particular path in 
the temperature-density plane crossing the glass transi- 
tion line: T e fj depends only on the final density. In this 
Letter we demonstrate that the findings of Di Leonardo 
et al. hold only for strongly correlating liquids (defined 
below) . We further argue and demonstrate that - for this 
class of liquids - from a single quench simulation one can 
predict the effective temperatures for any off-equilibrium 
jump. 

Reference documented the existence of a large class 
of liquids characterized by strong correlations between 
virial (W = PV — NksT) and potential energy (U) ther- 
mal equilibrium fluctuations at fixed volume, AW(t) = 
jAU(t). Strongly correlating liquids have a hidden (ap- 
proximate) scale invariance, which implies that they in- 
herit many - but not all - of the scaling properties 
of liquids interacting via inverse power-law potentials. 
Strongly correlating liquids include van der Waals liquids 
like Lennard- Jones type liquids but not, e.g., hydrogen- 
bonding liquids. Strongly correlating liquids have curves 
in their phase diagrams - "isomorphs" - along which sev- 
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eral static and dynamic properties are invariantP. These 
invariants derive from the fact that two microscopic con- 
figurations of two isomorphic state points, which trivially 
scale into one another, to a good approximation have 
proportional canonical probabilities. If the density is de- 
noted by p and the temperature by T, an isomorph is 
given by p 1 /T = Const. The exponent 7 - which may be 
slightly state-point dependent - can be calculated from 
a quench simulation utilizing the relation between the 
relaxing averages, (W(t)} = j(U(t)} + W . 

Because the canonical probabilities of scaled config- 
urations belonging to the same isomorph are identical, 
a jump between two isomorphic state points takes the 
system instantaneously to equilibrium (property (i)j^. 
Moreover, jumps from isomorphic state points to the 
same final state point show identical aging behavior 
(property In view of these properties the results 

of Di Leonardo et al. 5 may be understood as follows. 
A crunch from density pi to density p2 can be ideally 
decomposed into two parts (Fig. [I]): First the system 
jumps instantaneously from its initial state to the corre- 
sponding isomorphic state at the final density (i.e., the 
state which has the same p 1 /T as the initial state), see 
the blue curve of Fig. [l] This is an equilibrium state 8 . 
Thereafter the system begins to thermalize to the final 
temperature at constant density. If the crunch is made to 
a state with very high density, the thermalization takes 
extremely long or maybe infinite time, and the effective 
temperature may be determined from the FDT violation 
as detailed below. In this way any crunch corresponds 
to a quench to the final density with the same relaxation 
pattern. In particular, these two transformations should 
have identical FDT violation factors and identical effec- 
tive temperatures. 

These ideas should apply to any strongly correlating 
liquid, not just the single-component LJ system. To 
confirm this we simulated the standard Kob-Andersen 
binary Lennard- Jones (KABLJ) system™^. Following 
Di Leonardo et alW we subjected the KABLJ liquid 
to a number of instantaneous quenches and crunches 
and calculated the effective temperatures from the 
fluctuation- dissipation relation (FDR). Recall that for 
off-equilibrium systems the FDR in fc# = 1 units i^HHUl 
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FIG. 1: Patterns followed by the KABLJ liquid in different 
off-equilibrium density and/or temperature jumps. Consider 
for example the case of a crunch (horizontal dotted line), 
where the system is densified at constant temperature. This 
transformation is equivalent to a quench (right-most vertical 
line) from an isomorphic state point having the final density 
of the crunch. Thus the T e ff (filled circle) is identical for 
these transformations. In all processes represented here the 
liquid undergoes a glass transition characterized by an effec- 
tive temperature that can be measured from the fluctuation- 
dissipation relation (FDR) Eq. T e ff versus (final) density 
constitutes an isomorph, as discussed later in the text. 
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FIG. 2: Response vs. correlation function for several den- 
sity/temperature jumps for the KABLJ system. All FD-plots 
have fixed t = 10 4 (Monte Carlo steps) and t' varying from 
10 3 to 10 4 . All functions plotted here have the same reduced 
k- vector 10 (referring to the final density) (o) and the same re- 
duced microscopic In the crunch (o) we set |k| = 7.81 
corresponding to the reduced k- vector |k| = 6.78 (see^for de- 
tails). The crunch (o) overlaps very well with the quench (□) 
that takes the system from an initial state isomorphic to the 
one of the crunch to the same final state. Note also the good 
superposition of the additional quench (y) that takes the sys- 
tem from a state isomorphic to the initial state of (o) to a 
state isomorphic to the final one of (o). The full lines have 
the slopes predicted from the density scaling relation Eq. |2| 
for T e ff, adjusted only by a vertical shift to fit the data. 



Td t , X {t,t') = -X{t,t')d t ,C{t,t'). (1) 

Here C = (A(t)B(t')) is the equilibrium correlation func- 
tion of the variables A and B in the unperturbed situa- 
tion, the perturbing contribution to the Hamiltonian is 
SH = —eB, x(£, £') = (A(t)) / e\ e -to is the response of A 
to an infinitesimal external field of magnitude e applied 
at time t' < t, and X is the FDT violation factor. This 
is unity at short times (t — t f )/t' <C 1 while X < 1 in 
the long-time limit (t — t')/t r ^> 1. We chose as dy- 
namic variables A\^(t) = iV _1 J2j Vj exp(ik • rj(t)) and 
Bk(t) = N[Ak(t) + A_k(t)], where the sum is extended 
to all N particles of the system and rjj = ±1 is a random 
variable with zero mean. With this choice the correla- 
tion function studied is the self-intermediate scattering 
function, and C{t,t f ) is dimensionless. 

For quenches to low enough temperatures, at long 
times an effective temperature of the slow degrees of 
freedom T e ff is associated with the FDT violation fac- 
tor: T eff = T/J&EE®. The effective temperature reflects 
the slow structural rearrangements in the sense that the 
aging system behaves as if it were thermalized at T e f^. 
We obtained X by calculating the correlation function 
and the response function in the non-equilibrium regime 



by means of X = X(t) = -d[T X (t, t'))/dC{t, t')\ u which 
applies at long times (note that the correct X is found 
by taking this derivative at fixed t, not at fixed t^). 

Recently Berthier introduced a new method for cal- 
culating the response without applying any external 
field to an off-equilibrium Monte-Carlo simulation of the 
KABLJ^. We use his procedure and study the same 
dynamic variables for building the correlation function 
and the response as in Ref. |T7j Figure [2] shows the FD- 
plots for the KABLJ system during a number of tem- 
perature/density jumps. In Fig. [2] we test the construc- 
tion of equivalent crunches and quenches argued above: a 
crunch and a quench from initial isomorphic states (i.e., 
with the same p 1 /T) to the same final T and p (red cir- 
cles and squares). Clearly, the crunch overlaps well with 
the quench since - as argued above - they follow the same 
aging pattern. The exponent 7 was estimated by a linear 
fit of the parameter plot (W(t)) vs (U(t)) when the sys- 
tem is relaxing after a temperature jump from T = 2.55 
to T = 0.319 at fixed p = 1.264. The resulting value 
7 = 5.01 is used throughout this paper. Fig. [3] shows the 
linear relation that connects (W(t)) and (U(t)) when the 
system performs out-of-equilibrium jumps. 

Identical responses and correlations do not only ap- 
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FIG. 3: Average virial versus potential energy per particle 
during aging in two quenches for the KABLJ system. These 
quenches are performed between states with the same initial 
and final p 1 /T. At each time the off-equilibrium states are 
connected by an isomorph (e.g., the dashed line). In these 
two jumps the FDT violation factor X should be identical. 
The slope of this (W(t)) vs. (U(t)) plot is 7. The analytical 
equations used for drawing the isomorphs in the W-U plot 
here are reported in Ref.= 
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FIG. 4: Effective temperature as a function of density in 
several crunch and/or quenches (in double log-scale) for the 
KABLJ system. The effective temperature here is computed 
from the violation factor X, which is estimated as the slope 
that fits best the FD-plot (see Fig. [2J at points having 
(t — t')/t' > 1. The scaling exponent 7 is computed from 
potential energy-virial relaxation (see Fig. [3| as described in 
the text. The full line is the prediction of the density-scaling 
equation |2]) for T e ff. 



pear when a strongly correlating liquid is taken from two 
isomorphic state to the same state point. Supplement- 
ing properties (i) and (ii), strongly correlating liquids 
have a third interesting aging property (Hi): For two 
jumps (Ti,pi) (T 2 ,p 2 ) and (T 3 ,p 3 ) (T 4 ,p 4 ) be- 
tween mutually isomorphic initial and final states (i.e., 
p{/Ti = P3/T3 and p 2 /T 2 = P4/T4) the systems follow 
the same path in configuration space in reduced units 
(see^), and the dynamical equations governing the evo- 
lution of the particle trajectories are identical when cast 
in reduced units. Accordingly, the responses and corre- 
lations of two such jumps must be identical in reduced 
units. In Fig. [2] we show the reduced-unit C and \ °f a 
quench between initial and final states that are isomor- 
phic, respectively, to the initial and final states of the 
crunch described above (red triangles). The overlap be- 
tween the functions is good. This concept is visualized in 
Fig. [3j which shows the variables (W(t)) vs (U (t)) in two 
jumps (with the same initial and final p 1 /T) connected 
by an isomorph at each time t during the relaxation. 

A further consequence of property (Hi) is the follow- 
ing. Because the reduced-unit evolution of the dynam- 
ical properties is the same for the system in the two 
jumps, their FDR violation factors must also be identical, 
X2 = X4. Combining this equation with p 2 /^2 = P4/T4 
and expressing X via the effective temperature, we find 
p 2 /T e ff ?2 = pl/Tes^. This implies that 



p7T eff = Const. (2) 

This equation identifies the dynamic glass transition 
curve in the (T, p) plane defined in terms of the FDR 
effective temperature with an isomorph. This is consis- 
tent with the findings of Ref. [5] and the standard way of 
defining the glass transition, because the standard glass 
line in the (T, p) plane is located where the equilibrium 
relaxation time reaches a certain (high) value of order the 
inverse cooling rate. For strongly correlating liquids an 
isomorph is also an "isochronal" curve along which the 
(reduced) relaxation time is constantP. Figure [2] shows 
the slopes predicted by Eq. Q (lines); clearly the pre- 
diction is fulfilled. 

It is well known (see for example^ ancPS) that the ef- 
fective temperature is independent of the initial and fi- 
nal temperature if the initial temperature is high (the 
system is in a warm liquid state) and if the quench- 
ing temperature is low enough (i.e., in the linear regime 
where X = T/T e ff with constant T e ff). Equation Q pre- 
dicts the effective temperatures for all state points of the 
phase diagrams, with the exponent 7 as well as the con- 
stant found from one single aging experiment. In Fig. 
[4] we compare the T e ff measured in many crunches and 
quenches (not only involving isomorphic initial and final 
state points) with the prediction of Eq. Q . The agree- 
ment is very good. 

The above discussed simple aging properties are only 



4 



LENNARD-JONES GAUSSIAN MONATOMIC SYSTEM 




0.4 0.5 0.6 0.7 0.8 0.9 1.0 



C(t.t') 



FIG. 5: FD-plot in two off-equilibrium transformations for 
the MLJG system at the same final density and temperature. 
The potential defining this model is shown in the inset. In 
this case the two effective temperatures are very different. 



connected by the density-scaling relation. Moreover, for 
this class of liquids the effective temperature satisfies the 
density-scaling equation Eq. (|2|. Since the exponent 7 
and the constant of Eq. ([2| may be identified from a 
single quench simulation, the implication is that for a 
strongly correlating liquid the effective temperature of 
an arbitrary glass may be calculated from the results of 
a single jump simulation. 

The center for viscous liquid dynamics Glass and Time 
is sponsored by the Danish National Research Founda- 
tion (DNRF). 

I. AUXILIARY MATERIAL 

A. Simulation Details 

We simulate a binary mixture (80:20) of 1000 parti- 
cles interacting via Kob- Andersen Lennard- Jones poten- 
tial. Simulations are performed in the canonical (NVT) 
ensemble using the Metropolis Monte Carlo algorithm. 
Simulation parameters are taken from Ref. [9j 



expected to apply for liquids with isomorphs, i.e., 
strongly correlating liquids. To validate this we simu- 
lated the non-equilibrium dynamics of the monatomic 
Lennard- J ones Gaussian (MLJG) modeP^. The pair- 
potential of the MLJG has an additional Gaussian at- 
tractive well compared to the LJ liquid (see the inset 
of Fig. [5J; details about the model's potential and its 
glassy behavior can be found in Ref. 20, The MLJG 
liquid has WU fluctuations which correlate less than 2% 
at the state points studied here. As is clear from Fig. [5] 
two jumps to the same final density lead to quite differ- 
ent effective temperatures. Thus this system provides a 
counterexample to the observation by Di Leonardo et al. 
that the effective temperature depends only on the final 
density. 

We also investigated the relationship between the in- 
herent state energies in aging and at the equilibrium for 
the KABLJ and the MLJG (see 10 for more details). We 
find that only for the strongly correlating liquid KABLJ 
can one interprete T e ff as an indicator of which part of 
the energy landscape is visted during aging (as suggested 
by Sciortino et al. irP^). 

In conclusion, the existence of isomorphs for strongly 
correlating liquids explains the previously reported result 
for the LJ liquid that the effective temperature depends 
only on the final density of any jump (when tempera- 
ture and density are the externally controlled variables)^. 
We presented simulations of the aging dynamics of an- 
other strongly correlating liquid (the KABLJ system), 
as well as simulations of aging of a liquid without strong 
virial / potential energy correlations (the MLJG liquid). 
For a strongly correlating liquid it is always possible to 
produce equivalent density /temperature transformations 



B. Reduced Units in The MC Simulation 

The reduced length 8 is given by [=p 1 / 3 /. This defines 
the reduced k- vector k = p _1//3 k implying that, given a 
wave- vector ki at density pi , the k- vector with equivalent 
k at density p2 is k 2 = (p2/pi) 1 ^ 3 ki. 

The microscopic time r in the MC dynamics is set 
by the maximum random step attempted by the particle 
^max, i-e. t = r(S max ). This is fixed to minimize the 
relaxation time of the (equilibrium) correlation function 
at a given intermediate T and p (we set £ max at the same 
value chosen in I9|17l minimizing the decay-time of C at 
T = 0.75 and p = 1.2). It follows that two MC simu- 
lations (performed at pi and p2 respectively) have the 
same reduced time f if £ max is the same, that is when 
#max2 = (pi/p2) 1 ^ 3 ^maxi- However note that this is mi- 
nor adjustment of £ max : in the density range explored (in 
our simulations p varies from 1 to 1.53) t ax changes at 
most 15 %. 



C. Inherent States Energy: Equilibrium and 
Off-Equilibrium 

For a number of our simulations we evaluated the 
mean energy of the inherent states (local minima in the 
3N+1 dimensional potential energy landscape) visited at 
t = 10 4 . In Fig. [6] this is plotted against the mean inher- 
ent energy in equilibrium when the system is simulated 
at T = T e ff of corresponding aging experiment. For the 
strongly correlating KABLJ system there is good agree- 
ment between these two quantities, consistent with the 
interpretation that T e ^ indicates what part of the po- 
tential energy landscape is visited during the aging^. 
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KOB-ANDERSEN BYNARY 
LENNARD-JONES MIXTURE 
□ FromT=2toT=0.25 

(p=1.2039) T ef =0.855 
▼ From T=6.58 to T=0.835 

(,>=1 53) T eff =3.073 
• From T=2.55to T=0.22 

(p= 1.2039) T eff = 1.240 
A From T=2, p=1.2039 to 

T=0,25, p=1 .2642 T =1.733 
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INHERENT STATE ENERGY (in aging at f-10 4 ) 
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In contrast, this interpretation is clearly not valid for 
the non-strongly correlating MLJG system, supporting 
the notion that strongly correlating liquids have simpler 
physics than liquids in generaP. 



FIG. 6: Comparison between the (average) equilibrium in- 
herent states energy (at T = T e ff (£)) and the off-equilibrium 
inherent states energy (at time t) in several jumps in the glass 
(see legend) . For the KABL J systems (red symbols) these two 
quantities coincide, indicating that the aging system at t visits 
the same inherent states of the system in equilibrium at the 
temperature T = T e s(t) given by the FDR. This is clearly not 
the case for the MLJG systems (green symbols). The energy 
values of the MLJG are scaled by a factor 2.5 for an easier 
comparison with the KABL J energies. 
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